Disorder effects and current percolation in FeAs based superconductors 
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We report the influence of atomic disorder introduced by sequential neutron irradiation on the 
basic superconducting properties, flux pinning and grain connectivity. Two different polycrystalline 
Sm-1111 samples (SmFeAsOi-xFa;) and two Ba-122 single crystals (BaFe1.sCoo.2As2) were investi- 
gated. The monotonous decrease of the transition temperature with neutron fluence degrades the 
upper critical field, at least in the investigated temperature region. Pinning on the other hand is 
largely improved, with a different optimal defect concentration (fluence) in the two materials. The 
analysis of the current flow in the polycrystalline samples reveals weak link behaviour in the major- 
ity of grain connections and the existence of stronger grain connections. The density of the latter 
seems to be close to the percolation threshold (i.e. the minimum fraction for a continuous current 
path). Both types of connections are sensitive to disorder and degrade upon neutron irradiation. 
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I. INTRODUCTION 

The optimization of the superconducting properties, 
flux pinning, and intergranular current flow are impor- 
tant issues for applications of superconductors, which can 
be addressed by neutron irradiation studies^— . Fast neu- 
trons introduce atomic disorder and in many cases pin- 
ning efficient defects. The enhanced impurity scattering 
potentially changes the transition temperature 3 - - — , the su- 
perconducting energy gap^, the upper critical field 3 -^, 
and the magnetic penetration depth 9 . These defects are 
a model system for naturally grown or artificially pro- 
duced defects**^ and help to understand sample to sam- 
ple variations and to identify possible routes for material 
optimization 11 . Theoretical predictions can be checked^ 
and limitations of the current flow identifie d 11 ' 12 , since 
inter- and intragranular defects are differently affected. 
The influence of nanosized defects as bulk pinning cen- 
ters was unambiguously demonstrated^! 7 -. 

Fe-based superconductors offer promising supercon- 
ducting properties, i.e. comparatively high transition 
temperatures^— and high upper critical fields^ 7 - - — . 
The intergranular current flow, however, seems to be 
problematic^— . Very few irradiation studies on bulk 
samples^ ' 26 ' 27 and single crystals^ were performed until 
now. 

In this paper we report the influence of neutron irra- 
diation on the transition temperature, the upper critical 
field, and flux pinning. The inter- and intragranular cur- 
rent flow is investigated. 



II. EXPERIMENTAL 

One SmFeAsOi_ x F x sample (denoted as Sm-llllA in 
the following) was prepared at the National High Mag- 
netic Field Laboratory. The starting materials of As, Sm, 
Fe, Fe2C>3 and SmF3 were mixed and pressed into a pel- 
let, wrapped with Nb foil, and sealed in a stainless steel 



tube. The sealed sample was heat treated at 1160 °C 
for 6 hours in a high temperature isostatic press under 
a pressure of 280 MPa. The main phase of the sample is 
SmFeAsOi_ x F x , with a grain size of 10 to 15 microm- 
eters. The impurity phases include SmAs, SmOF and 
FeAs. The size of the sample was 1.1 x 1.7 x 3.7mm 3 . 

The second polycrystalline sample of SmFeAsOi_ x F x 
(sample Sm-llllB) was synthesized at the ETH 
Zurich^. SmAs, FeAs, SmF 3 , Fe 2 3 , and Fe were pul- 
verized and sealed in a BN crucible. At a pressure of 
3GPa, the temperature was increased to 1350-1450 °C 
within 1 h. After keeping this temperature for 4.5 h, 
the crucible was quickly cooled down to room tempera- 
ture. The sample investigated in this study was a piece of 
the pellet with an approximately triangular cross section 
(~ 2.5 mm 2 ) and a height of 1.9 mm. 

BaFe1.gCoo.2As2 single crystals were prepared by the 
self- flux method 30 under a pressure of 280 MPa at the 
National High Magnetic Field Laboratory. One crystal 
was characterized resistively, another inductively. Typi- 
cal dimensions of these crystals were 1.4 x 0.7 x 0.1mm 3 . 

Neutron irradiation was performed in the central ir- 
radiation facility of the TRIGA-Mark-II reactor in Vi- 
enna. The samples were sealed into a quartz tube and 
sequentially irradiated, starting with a fast neutron flu- 
ence (E > 0.1 MeV) of 4 x 10 21 m~ 2 . The highest cumu- 
lative fluence was 1.8 x 10 22 m~ 2 . 

Neutrons transfer their energy to the lattice atoms by 
direct collisions. The transferred energy must exceed the 
binding energy of the lattice atom to displace it, thus 
only fast neutrons lead to defects. The size of defects is 
expected to range from point defects to several nm, as 
observed in the cuprates 3 -^. Self shielding effects can be 
neglected, since the penetration depth of fast neutrons is 
of the order of a few centimeters, which is much larger 
than the sample dimensions. Only neutrons of low or 
intermediate energies are shielded efficiently because of 
large neutron cross sections at these energies. Some of 
these neutron capture reactions are followed by 7- or /3- 
cmissions, which may produce single displaced atoms 2 ®. 



These reactions are restricted to the surface in Sm-1111, 
since Sm is a paricular strong neutron absorber. The 
bulk will be penetrated only by fast neutrons. In con- 
trast, the recoil from 7- and /3-emissions might also con- 
tribute to the total defect density in the bulk of Ba-122. 
However, since we do not find any significant difference 
in the change of the transition temperature in the two 
compounds after irradiation, defects seem to result pre- 
dominantly from fast neutrons. 

The resistivities of sample Sm-llllA and of one Ba- 
122 crystal were measured at various fixed fields while 
cooling at a rate of lOK/h with an applied current of 
10 mA (single crystal: 300 /xA). Current and voltage con- 
tacts were made by silver paste (single crystal: silver 
epoxy). The distance between the two voltage contacts 
was around 1 mm for the polycrystalline samples and 
around 0.3 mm for the crystal. The transition temper- 
ature and the upper critical field were determined by 
means of a 90% criterion. 

Magnetization loops at various temperatures were 
recorded in a commercial 7T SQUID magnetometer 
(sample Sm-llllA) and a commercial vibrating sample 
magnetometer (Ba-122). The critical current density was 
calculated from the magnetization loops using the Bean 
model. A self field correction was made for the single 
crystal. The self field was calculated numerically, aver- 
aged over the sample volume, and added to the applied 
field /iqH, leading to J c (£)22. J c of sample Sm-llllA 
represents a rough estimation because of the uncertainty 
in the geometry of the current loops (grain size). The 
self field correction was therefore abandoned (J c (hqH)). 

The ac susceptibility at 33 Hz was measured in the 
same SQUID with an amplitude of 30 /xT (Ba-122 sin- 
gle crystal: 100 fiT) to determine T c and the shielding 
fraction. The demagnetization factor, which is needed 
for estimating the shielding fraction, was calculated nu- 
merically for the actual sample geometry. 

In addition, the remnant magnetic moment was 
measured (SQUID) as function of the maximally ap- 
plied field, -ffmax, in order to check for magnetic 
granularity^!" 3 ^!. 



III. RESULTS AND DISCUSSION 

Neutron irradiation introduces scattering centers, 
which increase the normal state resistivity. The normal- 
ized resistivity p n orm := p(55 K) / (p(300 K) — p(55 K)) in- 
creased from 0.49 in the unirradiated sample Sm-llllA 
to 0.55, 0.61, and 0.72 after irradiation to a fluence, F, 
of 4, 8, and 18 x 10 21 m~ 2 , respectively. The normaliza- 
tion to the phonon contribution (p(300 K) — p(55 K) = 
450 ±150 /ifkm) potentially cancels a reduction of the 
effective current carrying cross section due to secondary 
phases, voids, or cracks? 5 and errors in the determi- 
nation of the distance between the voltage taps. The 
nearly linear increase of p n orm with neutron fluence (~ 
1.3 x 10 _23 m 2 ) indicates that the defect density scales 
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FIG. 1. The transition temperatures of three samples as a 
function of the fast neutron fluence, F. The introduced disor- 
der decreases T c of all samples with a similar slope dT c /dF. 



with neutron fluence and that the resistivity is still far 
from saturation, where the mean free path of the charge 
carriers approaches the lattice parameter. The increase 
of impurity scattering is also evidenced by the resistiv- 
ity ratio p(300K)/p(55K), which decreases from 3.03 to 
2.83, 2.64, and 2.38. 

This resistivity ratio was only 2.2 in the Ba-122 sin- 
gle crystal and decreased to 1.9 after irradiation to 
4 x 10 21 m -2 . A clear increase of resistivity after irra- 
diation is thus observed in both compounds, in single 
and polycrystals. 

Neutron irradiation also decreases the transition tem- 
perature (Fig. [I} with a similar slope in both materi- 
als: dT c /dF = -9.1, -7.75, and -7.5 x lO^Km 2 
in Sm-llllA (resistive), Sm-llllB (inductive), and Ba- 
122 (inductive), respectively. This decrease of the tran- 
sition temperature caused by impurity scattering is a 
general property of superconductors with anisotropic 
energy gap 3 -^, and particularly with s± (extended s- 
wave 3 -^) gap symmetry 3 ^, or a consequence of multiband 
superconductivity!. 

Karkin et al2£ found a complete suppression of super- 
conductivity in La-1111 after irradiation to a fast neu- 
tron fluence of 1.6 x 10 23 m~ 2 , which is about one or- 
der of magnitude higher than the highest fluence of the 
present study. The slope dT c /dF found in our samples 
would predict a decrease in T c of between 12 K and 15 K 
at their neutron fluence. Either the decrease of T c with 
neutron fluence does not stay linear at high fluences, or, 
their sample behaves differently. Note that T c of the La- 
1111 sample was significantly more below optimum than 
that of our Sm-1111 samples. 

The upper critical field of sample Sm-llllA near the 
transition temperature is presented in Fig. [5] The de- 
crease in transition temperature shifts the B C 2(T) curves 
to lower temperatures, but the slope increases after irra- 
diation (see the insert). The interplay between the de- 
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FIG. 2. Temperature dependence of the upper critical field, 
B C 2, in sample Sm-llllA. The insert contain the same data, 
but the temperature was normalized by the transition temper- 
ature, T c . The influence of neutron irradiation is dominated 
by the decrease in transition temperature. Nevertheless, the 
slope increases after neutron irradiation, having its maximum 
after the first irradiation step (4 x 10 21 m~ 2 ). 



crease in transition temperature and the increase in slope 
suggests the presence of a maximum in B C 2 (0) at around 
4 x 10 21 m- 2 . 

In the Ba-122 system, the slope of B C 2(T) hardly 
changes after irradiation for H \\ ab leading to a general 
decrease of the upper critical field in this field orienta- 
tion. B c2 (T) becomes steeper for H || c, which results 
in a reduction of the upper critical field anisotropy (e.g. 
from 2.8 to 2.5 at 22 K)2&. Good agreement of the angu- 
lar dependence of B C 2 with anisotropic Ginzburg-Landau 
theory is found near Tc before and after irradiation. 

A. Flux pinning 

The critical current densities of the polycrystalline 
sample Sm-llllA and the single crystal Ba-122 are com- 
pared before and after neutron irradiaton to various flu- 
ences in Fig. [3] The single crystal data refer to a well 
defined crystallographic orientation of the applied field 
(H || c) and the induced currents (parallel to the ab 
planes), whereas the magnetic moment of the polycrys- 
talline sample results from currents flowing in all crys- 
tallographic directions providing only a rough estimate 
of the intragranular currents (see below). However, the 
magnitude and field dependence of J c in both samples is 
similar at 5K (upper panel), in particular after irradia- 
tion. Both materials show traces of a second peak (e.g. 
Ba-122 at 5K and Sm-1111 at 20 K), which disappear 
after irradiationS^. The second peak effect is generally 
assumed to indicate the transition from a low J c ordered 
vortex phase at low magnetic fields to a high J c disor- 
dered phase at higher fields. Introducing more disorder 
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FIG. 3. The critical current density increases in Sm-1111 
(sample A, solid symbols) and Ba-122 (open symbols) due to 
the introduction of pinning efficient defects by neutron irra- 
diation. The x-axis refer to B for the Ba-122 crystal and to 
fioH for the Sm-1111 sample, since no self field correction was 
made in the latter case. J c is an estimate for the intragranular 
currents in Sm-1111. 

shifts the transition field to lower values and eventually 
leads to the disappearance of the ordered phase, as is ob- 
served in our samples upon irradiation. The introduced 
defects are efficient pinning centers in both materials. 

The self field J c increases after the first irradiation 
step, but does not change much upon further irradia- 
tion. Higher fluences decrease the field dependence of J c 
in Sm-1111, but not in the Ba-122 single crystal, where 
the second irradiation step decreases J c at 20 K, a conse- 
quence of the reduced transition temperature. 

The temperature dependence of the critical current 
densities at 1 T is shown in Fig. |4j The relative en- 
hancement after irradiation (by up to one order of mag- 
nitude) increases with temperature (sample Sm-llllA, 
solid symbols), which indicates that the pinning energy 
of the radiation induced defects can better compete with 
the thermal energy than the pinning energy of the as 
grown defect structure. A reasonable explanation is a 
larger defect size, since J c is not very sensitive to the 
defect density (fluence). 

In the Ba-122 crystal on the other hand, the temper- 
ature dependence hardly changes after irradiation, only 
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FIG. 4. Critical current densities at IT in samples Ba-122 
and Sm-llllA as a function of temperature prior to and after 
neutron irradiation to various fluences. 











/ 

□ 








■ 


5 K 


35 K 




/ 


■ 


— □— 


Sm-1111A ■ 


/ 

,80888888=8=8=0= 


— •— 


— o— 


Sm-1111B 

^*=8=8=4 







0.00 0.02 0.04 0.06 0.08 0.10 



S(T) 



FIG. 5. Ac susceptibility as a function of the applied dc field. 
While the global currents are rapidly suppressed in sample 
Sm-llllA, they remain significant (>• 2.5 x 10 4 Am~ 2 ) in 
sample Sm-llllB. 

at high temperatures the slope becomes slightly steeper, 
which is caused by the reduced transition temperature 
(20 K is close to T c ) . The pinning energies of the as-grown 
and the radiation induced defects seem to be similar. 



B. Intra- and intergranular currents 

The global current flow was checked by ac susceptibil- 
ity measurements (Fig. [5]). The ac susceptibility, x'> °f 
sample Sm-llllA rapidly increases (shielding decreases) 
with applied dc field, which indicates a rapid decoupling 
of the grains. A metallic Fe-As wetting-phase, which 
can be passed by the supercurrents only via Josephson 
tunneling, was found at the grain boundaries of a compa- 



rable sample 2 ^.. These superconducting tunnel currents 
are known to be rapidly suppressed by magnetic fields. 

Shielding is perfect = — 1) within experimental ac- 
curacy in sample Sm-llllB. This means that the super- 
currents only flow at the sample surface and the pen- 
etration depth of the applied ac field (6 = fP c /J c )is 
negligible compared to the smallest sample dimension 
R s ; thus, the intergranular currents must be much larger 
than H ac /R s « 2.5 x 10 4 Am" 2 . The field independent 
susceptibility proves the existence of significant global 
currents, at least in the field region relevant for the mea- 
surement of the remnant magnetic moment. This does 
not necessarily mean that weak (Josephson) grain cou- 
pling is absent in this sample. However, the number of 
stronger grain connections must be significantly above 
the percolation threshold 39 , which is defined as the min- 
imum fraction for a continuous current path. 

In order to investigate the influence of disorder on the 
inter- and intragranular currents, the remnant magnetic 
moment, m rom , was measured as a function of the max- 
imally applied field, H max (e.g. solid circles in Fig. [Bb). 
The magnetic moment was measured at each field i? max 
(Fig. [BJi) and after reducing the field to zero (circles in 
Fig.EJj). If all current loops, which were induced by the 
application and removal of the external magnetic field, 
are of similar geometry, a peak in the derivative of the 
remnant magnetic moment with respect to the logarithm 
of the maximally applied magnetic field, 

d^rem /-.\ 

dlog(ff max )' ( ' 

occurs, when ff max approximately equals the Bean pen- 
etration field H*—. The critical current density can be 
estimated from H* w J C R, if the representative geometry 
of the current loops, R, is known (typically the smallest 
sample or grain dimension). 

Two peaks are observed in granular materials^, if 

jintra^g < jinter RS ) where jintra and jintcr denote the 

critical current density within or between the grains, re- 
spectively, and R s (R £ ) is the smallest sample (grain) 
geometry. Flux first penetrates the whole sample along 
the grain boundaries, while the inner parts of the grains 
are still free of vortices. If the above condition is not 
fulfilled, the peaks overlap and cannot be properly dis- 
tinguished. 

The logarithmic derivative ([T]) of sample Sm-llllB is 
presented in Fig. [7J Two peaks occur at 20 and 35 K (also 
at 12.5 and 27.5 K, not shown) in the unirradiated state 
(solid squares). The low field peak is usually ascribed to 
intergranular currents shielding the whole sample. This 
peak is located at about 14 mT at 5K (solid squares in 
Fig. Ek), which corresponds to a critical current density 
of around 10 7 Am -2 decreasing to about 5 x 10 6 Am -2 
at 35 K. 

It is a priori not clear whether the low field peak sim- 
ply arises from a rapid flux penetration immediately after 
ff max exceeds the lower critical field, of the sam- 

ple, since the reversible magnetization strongly changes 
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FIG. 6. Magnetic moment as a function of the (increasing) ap- 
plied magnetic field A/io-ffmax (panel a). Measurements were 
done in a background field, HoHo, of 0T (circles) and 0.1 T 
(triangles). The field was reduced to fj,oHo between each value 
of /io^max in order to measure the remnant magnetic moment 
jn rcm (panel b). A peak occurs in dm rem /d(log Ai/ max ), each 
time Af/max approximately equals a penetration field (panel 
c). The low and high field peaks result from flux penetraion 
into the whole sample (intergranular peak) and the grains 
(intragranular peak), respectively. All peaks shift to lower 
fields at noHo = 0.1 T. The shift of the intergranular peaks 
is more pronounced, indicating a stronger field dependence of 
jmtcr pj a t a were obtained at 5 K on sample Sm-llllB after 
the first irradiation step (4 x 10 21 m -2 ). The inserts enlarge 
interesting details. 



around H c ±. In order to exclude this scenario, the mea- 
surement was repeated after field cooling the sample at 
0.1 T. The measurement sequence was repeated with this 
constant offset field, /j,qHo — 0.1 T. The field was in- 
creased by AH max and decreased to 0.1 T, where the 
trapped magnetic moment was measured. The results 
at 5 K are compared to the standard measurements 
(PoHq — T) in Fig. [5] (the moment measured imme- 
diately after field cooling was subtracted). The low field 
peak persists, but shifts to lower fields (by a factor of 
1.7), which is a result of the field dependence of the crit- 
ical currents. As can be seen in the insert of panel b, the 
remnant magnetic moment is larger (between Ai/ max = 2 
and 8 mT) first at 0.1 T, which results from an earlier flux 
penetration (cf. panel a). Then the larger critical cur- 
rents at T lead to a larger trapped flux, until the same 
behaviour is repeated near the high field peak (m rem at 
0.1 T first larger then smaller than in zero field). In any 
case, the low field peak is well above the first penetration 
field and obviously not caused by the reversible magneti- 
zation, which confirms its interpretation as an intergran- 
ular peak. 

Neutron irradiation shifts the intergranular peak to 
lower fields (Fig. [7J), which indicates a degradation of 
the intergranular currents. The effect is small after ir- 
radiation to 4 x 10 21 mT 2 (circles), but significant at 
1.4 x 10 22 m~ 2 (line graph). Such a degradation of the 
intergranular currents after neutron irradiation was also 
found in Bi-tape o 12 i 41 and coated conductors^ 2 -. 

The high field peak usually originates from the intra- 
granular currents shielding only individual grain o 20 i 34 ' 40 . 
In sample Sm-llllB, however, the high field peak splits 
at 5K, which sheds some doubts on this standard in- 
terpretation. A similar peak split at 5 K was reported 
for a polycrystalline Nd-1111 sample 2 ^. Three peaks are 
clearly visible at all temperatures after neutron irradia- 
tion (solid line in Fig. EJ). 

It will be argued in the following, that the intermediate 
field peak is caused by clusters of better connected grains 
and that the high field peak is the intragranular peak. 
These two peaks are only visible at 5 K in the unirradi- 
ated sample and merge at 12.5, 20, 27.5 and 35 K. This 
temperature dependence excludes the scenario, where the 
two peaks are caused by currents of the same magnitude 
flowing on different length scales (e.g., grains and clus- 
ters, or grains (clusters) of different sizes). This is con- 
firmed by the changes after irradiation, which are quite 
different for these peaks. While the high field peak shifts 
to higher fields, the intermediate peak shows a less sys- 
tematic behavior. It shifts to lower fields at 5K and 
12.5 K but approximately retains its position at higher 
temperatures (Fig. [7]). The shift of the high field peak 
obviously results from improved pinning (larger intra- 
granular currents, see Sec. MI AT) which also enhances the 
hysteresis (e.g., larger saturated remnant moment). The 
decrease of the intermediate peak field at low tempera- 
tures after irradiation indicates that the corresponding 
currents are limited by grain boundaries. This scenario 
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FIG. 7. Evolution of the peaks in dm rcm /d(log J/ max ) (cf. 
caption in Fig. [6] here with fioHo = T) with temperature 
and neutron fluence. The irradiation decreases the intergran- 
ular currents (low and intermediate field peak). The intra- 
granular (high field) peak shifts to higher fields, as a conse- 
quence of the increase in J^ ntra . The intermediate and high 
field peak merge at temperatures above about 10 K in the 
unirradiated sample (Sm-llllB). 

is supported by the field dependence of the peak posi- 
tions (Fig. While the high field peak shifts only by 
about 5%, the low and intermediate peak fields decrease 
by a factor of roughly two due to the increase of the 
background field from (solid circles) to 0.1 T (open tri- 
angles) . (The intermediate field peak degrades to a bump 
at0.1T.) 

Only a rough estimate can be given for J c in the current 
loops which are responsible for the intermediate peak. 
Since the size of the loops must be smaller than the sam- 



ple geometry and larger than the individual grains, the 
corresponding current density is between 2 x 10 8 Am -2 
and 2 x 10 9 Am' 2 at 5K and 0T. 

At least two types of grain connections exist in our 
polycrystalline samples: (i) the wetted grain boundaries 
(type- A), whose currents are rapidly suppressed by mag- 
netic fields and (ii) grain boundaries over which cur- 
rents can pass also at high magnetic fields (type-B). The 
number of type-B connections is below the percolation 
threshold in sample Sm-llllA and above the percola- 
tion threshold in sample Sm-llllB, as indicated by the 
ac susceptibility. The rapid decrease of x' directly shows 
the suppression of global currents after "switching-off" 
type-A connections by applying an external magnetic 
field. It was pointed out in^ 2 - that samples of appar- 
ently similar microstructure (probed by x-ray and SEM) 
showed a large variation in the intergrain current density. 
It is therefore reasonable to assume that the number of 
type-B connections is close to, but below the percolation 
threshold in our sample Sm-llllA. In fact the low field 
peak was not present in this sample, since the global cur- 
rents were too small, but we observed the intermediate 
and high field peaks (at comparable fields) which might 
just reflect the structure of the incipient spanning cluster 
at the percolation threshold. In the "nodes, links, and 
blobs" picture, the incipient spanning cluster is composed 
of blobs which are connected by links^. Removing one 
link decomposes the spanning cluster, blobs are multi- 
ply connected sites (clusters of type-B connected grains). 
The two peaks in sample Sm-llllA could correspond to 
flux penetration into the blobs and grains respectively. In 
a comparable sample which obviously contained a span- 
ning cluster of type-B connections, the links were made 
directly visible by scanning laser microscopy 2 ^. 

In sample Sm-llllB, the density of type-B connections 
seems to be significantly above the percolation threshold, 
since shielding of very small ac fields was nearly perfect 
(see above). However, the clusters of type-B connected 
grains may persist well above the percolation threshold 
and be responsible for the intermediate field peak. 

IV. CONCLUSIONS 

Atomic scale disorder is not a promising way for im- 
proving the basic superconducting properties (-B C 2, T c ) in 
the Fe-based superconductors. Bulk pinning on the other 
hand is strengthened by the introduction of nanoscale de- 
fects. The critical current densities in the grains gener- 
ally increase after the first irradiation step (fast neutron 
fluence: 4 x 10 21 m -2 ), but start to decrease again upon 
further irradiation in the Ba-122 single crystal, while a 
further improvement is found at high fields in Sm-1111 
up to a fluence of 1.8 x 10 22 m~ 2 . 

We found indications that the majority of grain con- 
nections are weak links which decouple in magnetic fields. 
We identify a second type of connections over which su- 
percurrents flow at high fields in a percolative manner, 
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since their density seems to be close to the percolation 
threshold. 
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